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ABSTRACT 

We use a new set of cold dark matter simulations of the local universe to investigate the dis- 
tribution of fossils of primordial dwarf galaxies within, and around the Milky Way. Throughout, 
we build upon previous results showing agreement between the observed stellar properties of 
a subset of the ultra-faint dwarfs and our simulated fossils. Here, we show that fossils of the 
first galaxies have galactocentric distributions and cumulative luminosity functions consistent 
with observations. In our model, we predict ~ 300 luminous satellites orbiting the Milky Way, 
50 — 70% of which are well preserved fossils. Within the Milky Way virial radius, the majority 
of these fossils have luminosities Ly < 10 5 L Q . Despite our multidimensional agreement with 
observations at low masses and luminosities the primordial model produces an overabundance 
of bright dwarf satellites (Ly > 10 4 L©) with respect to observations where observations are 
nearly complete. The "bright satellite problem" is most evident in the outer parts of the Milky 
Way. We estimate that, although relatively bright, the primordial stellar populations are very 
diffuse, producing a population with surface brightnesses below surveys detection limits and 
are easily stripped by tidal forces. Although we cannot yet present unmistakable evidence for 
the existence of the fossils of first galaxies in the Local Group, the results of our studies sug- 
gest observational strategies that may demonstrate their existence, i) The detection of "ghost 
halos" of primordial stars around isolated dwarfs would prove that stars formed in minihalos 
(M < 10 8 M Q ) before reionization, and strongly suggest that at least a fraction of the ultra- 
faint dwarfs are fossils of the first galaxies, ii) The existence of a yet unknown population 
of ~ 150 Milky Way ultra-faints with half-light radii r^i ~ 100 — 1000 pc and luminosities 
Ly < 10 4 L Q , detectable by future deep surveys. These undetected dwarfs would have the 
mass-to-light ratios, stellar velocity dispersions and metallicities predicted in this work. 



1. Introduction 



Over the last decade, since iKlypin et all (119991) and iMoore et al.l (119990 showed that the number of 
dark matter subhalos expected around a Milky Way mass halo is two orders of magnitude above the num- 
ber of known satellites, significant effort has been made in observation and theory to solve the substructure 
problem in cold dark matter (CDM) cosmology. Observational discoveries have redefined the 'Missing 
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Galactic Satellite Pr oblem' with the number of observed Milky Way satellites now somewhat closer to theo - 



The discovery of the u ltra-faint dwa rfs (Belokurov et al. 



retical expectations (iSimon & Gehal 120071: iTo llerud et al. 2008; Bovill & Ricotti 2009; Macci o et al 



Willman et al1l2005al lbl: Izucker et al]|2006al Jbl: iGeha et al. 



2006, 120071 : Ilrwin et all 120071 : Walsh et al 



01. 



201 



2007; 



2009) has roughly doubled the known Milky Way 



and M3 1 satellite populations since 2004. However, observations alone cannot fully explain the discrepancy 
between luminous satellites and CDM substructure. 

One of the core theoretical issues of the "missing satellites" problem remains the relationship between 
luminosity of satellites and the virial mass at formation of their dark matter halos. From a theoretical 
perspective, the fundamental question to b e answered is : "What is the minimum halo mass which can 



host a luminous galaxy?" Previous studies jEfstathiou 



Venkatesan et al 



2001; Ri cotti & Ostrik er 2004; 



1992 



Ricotti et al 



Tho ul & Weinberdl 19961 : Bullock et al.ll2001 



2005) have shown that this critical mass is 



set by the reheating history of the intergalactic medium (IGM) (|Ricotti et all 120001) and by feedback loops 
operating before reionization which determine whether low mass minihalos with M < 1O 8 M are able to 
accrete gas from the IGM and form stars. Simulations show that for halos with masses M < 10 s M the 
local and stochastic components of galaxy feedback produce mini halos with the sam e dark matter mass, but 
with stellar masses which vary by several orders of magnitude (IRicotti et all 120081) . It is likely unjustified 
to assume a sharp mass threshold separating dark and luminous halos and/or a tight relationship between 
dwarfs' luminosities and their total mass, at the faint end of the luminosity function. This is one of the main 
motivations for the present study. 

Cooling in halos with masses at formation, M > 10 8 M {v max > 20 km/s) is initiated via read- 
ily available hydrogen Lyman-a emission. However, cooling in minihalos with T v i r < 10 4 K requires 
the formation of either molecular hydrogen or pre-enrichment with metals from nearby galaxies. The bal- 
ance between the destruction and formation of H2 and metal transport in the IGM governs whether or not 
low mass minihalos can initiate cooling and form stars. If ultraviol et H2 dissociating radiation dominates , 
star formation in the f irst minihalos may be suppressed or delayed jHaiman et alJboOO ; Ciardi et al. 200C ; 
Machacek et aDboodlwise & Abelll2007l : lo'Shea & Normanll200^ : lMaio et alJboidl) . However, if ionizing 



ultraviolet radiation {hu > 13.6 eV) is the dominant feedback mec hanism, H? format i on can be catalyzed 
inside relic HII regions and on the edges of Stromegen spheres ( Ricotti et al. 2001 . 2002allbl: lAhn et al. 



2006: IWhalen et al.ll2008l : IWise & Abelll2008h . allowing star formation to be more widespread in minihalos 
with mass M 10 8 M , before reionization. It is important to emphasize that, regardless of the interplay 
between the creation and destruction of H2, not every minihalo will host a galaxy. 

The ultra-faint dwarfs are an excellent laboratory for testing models of star formation in the lowest 
mass halos, due to their extremely low luminosities: 10 2 — 10 5 L , low dynamical masses: M 10 7 M^ 



2010) , and their est imated large number ( Ko posov et al. 



dstrigari et alJl2008l:Jwalker et allboidlWolf et al 

2008 : lTollerud et alJboOslJwalsh et al-hood l. In lBovill & Ricottil (I2009T) and Bovill & Ricotti (2010a) (here- 



after Paper I), we argue th at the properties of the ultra-faint dwarfs are consistent with those of simulated 
primordial fossil galaxies. iKravtsovl (|20ld) has argued that the missing galactic satellites can be accounted 
for by a simple model in which the fraction of baryons turned into stars, /* decreases with the virial mass in 
a halo before its infall into a Milky Way. Not surprisingly, our simulations of pre-reinization dwarf galaxies 
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naturally prod uce such a decreas e of star formation efficiency with mass, although with increasing scatter at 
small masses (IRicotti et alj |2008). 



Assuming the o bserved ultra-faints belong to an isotrppical ly distributed satellite population around the 
Milky Way (but see, Metz et al. 2007 . 2009 : Bailin et aDboOsh . one infers from observations a minimum 
of 60-65 satellites within 200 kpc of the Galactic center. This number is consistent with all the satellites 
being subhalos that formed after reionization. There are about 90 (60) subhalos that during their evolution 
had a maxim um circular velocity a bove 20 km/s (30 km/s), thus likely formed most of their stars after 
reionization (IBovill & Rico tti 2009) (hereaft er BR09). However, w e know that most satellites with Ly < 
10 4 L Q are undetectable beyond 50—100 kpc ( Simon & Gehall2007 ). Once luminosity corrections are folded 
into estimates of the satellite counts, t he number of satellites inferred from observations is ~ 200 — 250 
(|Tollerud et allbood : Iwalsh et al.ll2009h . sugg esting that aro und the Milky Way there could be ~ 100 - 150 
true pre-reionization fossils (for a review see. iRicottil 120 1 Oh . Nevertheless, these luminosity corrections are 
uncertain and model dependent, hence it is largely unknown how many ultra-faints may exist in the outer 
parts of the Milky Way. Indeed, the dimmest ultra-faint dwarfs within 50 kpc of the Milky Way may belong 
to a different population: according to our simulations, they have properties that are inconsistent with well 
preserved fossils, suggesti ng that may have been shaped by tidal forces. Some observations also seem t o 
reinforce this hypothesis jMcGaugh & WolJ20ld : lwillman et aDboiolJSand et all2009l : lFrebel et allioioh . 
In summary, although there is circumstantial evidence that a fraction of the new ultra-faint dwarfs represent 
a well preserved primordial population in halos with infall masses corresponding to v max ^ 20 — 30 km s~ 2 , 
we have yet to find observational proof of their existence. 

But is this "primordial scenario" consistent with observations when we move beyond the virial radius of 
our own galaxy and peer into the voids? Unless the local component of galaxy feedback is very strong, star 
formation should proceed similarly in small mass halos regardless of a halo's location relative to the Local 
Group. Therefore, the voids should be populated with luminous objects (see Figure 1 of Paper I). As first 
noted in IPeeblesi (1200 ll) . they are not. The number of dwarf galaxies with absolute ma gnitude My > — 1 6 



(Ly < 2 x 1 8 L^) observed i n the voids is s maller than expected i n CDM cosmologies (IKarachentsev et al. 



2004, 120061 : iTully et all 120061) . According to ITikhonov & Klypinl (|2009r) . the luminosity function can only 
be reconciled if halos with v max < 35 km s _1 are dark. However, such a large mass threshold for star 
formation would produce less than 35 luminous satellites between 100-200 kpc of the Milky Way. That 
is inconsistent with observations unless satellites with Ly < 10 4 L do not exist beyond 100 kpc (see 
Figs. 3-4 in this paper). 

Alternatively, if the star formation rate is primarily determined by halo mass, the void phenomenon 
can be reconciled with CDM by using a halo occupation distribution for which the M/L ratio increases 
with decreasing halo m ass. However, this solution has only been tested for halos with M > 10 10 M Q 
(|Tinker & Conroyll2009r) . three orders of magnitude more massive than our fossils. When we extend their 
Mj L ~ M~ l relation to our mass range we obtain a Mj L ratio ^ 10 5 for the fossil population (M £ 10 s M Q ). 
This would produce an "ultra-faint" p opulation with M/L ratios ~ 2 — 3 orders of magnitude greater than 
those seen for the primordial fossils in lRicotti et al. ( 2002b ) 
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In this paper, we investigate the following conundrum: can we simultaneously account for the predicted 
and observed subhalo population around the Milky Way and the lack of isolated galaxies with My > —16 
in the voids? Is a primordial origin scenario for ultra-faint dwarfs, consistent with the number of dwarfs in 
the voids? The simulations described in Paper I allow us to address both of these populations within the 
same theoretical framework. As a result of this study, we will propose three new observational strategies to 
search for evidence of fossil galaxies in the local universe. We will argue that, if stars formed in halos with 
masses M < 1O 8 M there should be diffuse "ghost halos" of primordial stars around isolated dwarfs, in 
addition to the undetected population of lower surface brightness ultra-faint dwarfs. 

This paper is structured as follows. § [2] and [3] summarize the numerical method and treatment of the 
observations, covered in detail in Paper I. In §[4] we discuss the radial distribution of the fossils of the first 
galaxies and make comparisons to other N-body simulations (§ 14.11 ) and we show agreement between the 
radial distributions and cumulative luminosity functions of primordial fossils in our simulations and the 
observed ultra-faint satellites (5 I4.2I ). We also discuss a possible origin of Willman I type ultra-faint dwarfs 
and the exception of Segue I. In § |5l we describe an overabundance of bright satellites in the outer parts of 
the Milky Way halo and investigate possible solutions. In §|6l we propose observational tests for our models. 
Discussion and a summary are presented in § [7] 



2. Numerical Method 

In this paper, we continue our analysis of the simulations described in Paper I. In brief, we used a set 
of hybrid initial conditions for a ACDM N-body simulation to give us the resolution necessary to study 
the primordial fossils on Local Volume scales. Our high resolution region, centered on a 'Local Group,' is 
generated as follows. Instead of using a uniform grid of particles, on scales I < 1 M pc, the positions and 
velocities of our halos are set by the final outputs of the lRicotti et al. (l2002alJbLl2008h simulations, hereafter 



referred to as the pre-reionization simulations. Perturbations on scales / > 1 Mpc are added to the initial 
positions and velocities of particles at z = Zi n u using a coarse resolution simulation run from z = 40 to 
z = 0. For more detail, we direct the reader to the appendix of Paper I. In addition to the non-uniform 
initial grid of particles, the particles in the high resolution region have masses set by the mass function of 
the final pre-reionization outputs. Thus, we do not need to resolve these pre-reionization halos, just trace 
their merger history, tidal disruption, and positions from reionization to present day. 

Paper I details the two realizations of our initial conditions Our first order realization (see Appendix 
in Paper I for details) uses the same pre-reionization output across the entire high resolution region, and 
consequently does not account for the slower rate of structure evolution in the voids. In contrast, our second 
order initial conditions use four outputs from the pre-reionizat ion simulations and do account for slower 



rate of structure formation in the voids (|Barkana & Loebl 120041) . Near the Milky Way, where this paper is 



concerned, we find no substantive difference between the two realizations. We therefore focus our discussion 
on the second order initial conditions. 



All the simulations discussed in this paper were run using Gadget 2 (ISpringell 120051) and analyzed 
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using the halo finder AHF (IKnollmann & Knebd 120091) . Simulation and data analysis were run on the 
beowulf cluster Deepthought at the HPCC at the University of Maryland. At z = we define a bound 
halo as anything identified by AHF that, for the parameters we used, produces mass functions complete for 
halos with a number of particles N > 50. Paper I shows that our hybrid initial conditions and non-uniform 
particle masses produce results consistent with traditional CDM N-body simulations on both Local Volume 
and galactic subhalo scales. The parameters of the three 'Milky Ways' contained within our two highest 
resolution runs are listed in Table Q] 



2.1. Fossil Definition 

Dark matter halos identified in our simulations at z = are divided into three populations based on 
their ability to accrete gas from the IGM and form stars after reionization. A just virialized halo is able to 
accrete gas after reionization only if its maximum circular velocity, v m ax, is larger than a filtering velocity, 
Vfut- In this work, we use Vfm = 20 — 30 km s _1 , corresponding to the threshold for cooling via Lyman- 
a. Paper I and GK06 show that the exact choice of the filtering velocity does not significantly change the 
results. If a z = halo has v max > Vfm, we classify it as a non-fossil. In the modern epoch, non-fossils can 
be identified as dlrrs which have been accreting gas and forming stars continuously since reionization. 

Any halos whose present day v max is below the filtering velocity is a candidate fossil. A candidate 
fossil for which v max > Vfut at any point during its evolution may have accreted gas and/or formed stars 
after reionization. We include them in the "non-fossil" group, as post-reionization stars may be the domi- 
nant stellar population. When we need to distinguish these "candidate fossils" from the broader "non-fossil" 
group, we will refer to them as polluted fossils. Finally, any candidate fossil for which v max < v fm from 
reionization to the modern epoch is a true fossil. True fossils formed the majority (> 70% ) of their stars be 



fore r eionization and today would be relatively diffuse systems of old stars devoid of gas (IRicotti & Gnedin 



2005). The term fossil will only refer to "true fossils". Unless otherwise specified, the term non-fossil 
will apply to any halo which could have accreted gas after reionization regardless of its maximum circular 
velocity at z = 0. 

The division of the Milky Way and M31 satellites into fossils and non-fossils is shown in Table 2. We 

Table 1. Table of 'Milky Ways' 



Name Run Mass R v i r v max 
(1O 12 M ) (kpc) (km s" 1 ) 



MW.l 


C 


1.82 


248 


203 


MW.2 


D 


0.87 


222 


196 


MW.3 


D 


1.32 


194 


177 
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distinguish between the seven, classical RG05 fossils above the 10 Lq threshold, the RG05 fossils with 
Ly < 1O 6 L and the ultra-faint dwarfs discovered since 2005. 



3. A Note on Observations 



We approa ch the observ ations as follows. The majority of the information on the classical dwarfs 
comes from the iMateol (119981) review. For the ultra-faint dwarfs we g enerally defer to m easurements with 
the smallest error bars with some weight given to more recent work (IWalker et al.1 (120091) ). We direct the 
reader to Paper I of this series and BR09 for a more complete discussion of these criteria. 

When calculating the observed distributions of dwarfs arou nd the Milky Way , we account for two 



effects, the sky coverage of the SDSS, and its detection efficiencies (Walsh et al.112 009: Ko posov et al . 2008). 
For the classical dwarfs, we assume the entire sky has been covered and only apply sky coverage corrections 
to the ultra-faint population. To correct for the SDSS sky coverage, we assume that the satellite distribution 
around the Milky Way is isotropic, and multiply the number of ultra-faints by 3.54 to account fo r the nearly 
three-quarters of the sky not surveyed by SDSS, now past Data Release 7 (Abazaiian et alj|2009). However, 



brigh t satellites of the Milky Way are distribut ed very anisotropically (IKroupa et al.l 12005 : 



Zentner et al 



2005: LMetz et al.ll2007l l2009i : iBozek et al.ll201ll) . so the assumption of isotropy may not be a good one. 



Next, we apply a correction for the detection efficiency of the SDSS using the results from I Walsh et all 



(|2009l ). If an ultra-faint is bright enough to be detected with 99% efficiency, we assume the sample is 
complete for that luminosity and distance. However, if the ultra-faint is too dim for 99% detection, but 
bright enough to be detected half the time, we assume that, statistically, there is another satellite with similar 
luminosity and distance missed by SDSS. This second correction produces only a minor increase in the 
number of satellites; approximately one additional satellite over a total of ~ 60 from sky coverage correction 
alone. 

As discussed in Section 4 of Paper I, we also divide the ultra-faint dwarfs into two groups. The first 
is a group of seven, including CVnl and II, Hercules, Leo IV, Leo T, and Ursa Major I and II, have half 
light radii and surface brightnesses which are consistent with the stellar properties of fossils of the first 
galaxies. In contrast, the second group composed of five members, including Willman 1, Segue 1 and 2, Leo 
V, and Pisces II, have half light radii which are too small, and surface brightnesses which are too high, to be 
consistent with the simulated primordial population. 

Table [3] shows the Milky Way ultra-faint dwarfs classified in two ways, (i) the satellites within and 
beyond 50 kpc, and (ii) those which are consistent (opened circles in right panel of Figure 11 of BRlla) 
and those which an inconsistent (filled green circles in the same figure) with the expected properties of 
fossil dwarfs. Excepting Leo V and Pisces II, the dwarfs which have half-light radii significantly smaller 
than the simulations are also within 50 kpc of the Milky Way. Conversely, all the ultra-faint dwarfs which 
are consistent with the simulations are beyond 50 kpc, except Coma Ber.. In this paper, we will refer to 
this second population collectively as the "inner ultra-faints," to emphasize their location within the Local 
Group. 
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Fossils 






Non-fossils 


RG05 only 


RG05 &BRlla 


BRllaonly 




LMC 


Sculptor 


Draco 1 


Bootes I k II 




NGC55 




Phoenix 


CVn I k II 




Sextans A k B 




Sextans 


Hercules 




SMC 




Ursa Minor 


Leo IV k T 2 


Milky Way 


WLM 

Carina 
Fornax 
GR8 
Leo 1,11 & A 
Sagittarius 






Pisces II 




IC 10 


And I k H 


AndV 


And XI XII 




IC 1613 


And III 




And XIII k IV 




IC5152 


And VI 




And XV k XVI 




M32 


Antila 




AndXVII k XVIII 




NGC 185 


KKR25 




And XX k And XXI 


M31 


NGC 205 
NGC 3109 
NGC 6822 
DDO 210 
LGC3 
Pegasus 






And XXII k XXIII 
And XXIV k XXV 
And XXVI k XXVII 


Isolated 






Cetus 3 
Tucana 





Table 2: Table of all the Local Group dwarfs divided by their host (or lack there of) and their non-fossil or 
fossil status. We have further divided the fossils into three groups, those considered fossils in RG05, but 
with Ly > 1O 6 L , the "classical" fossil dwarfs with luminosities below the 1O 6 L threshold from BRlla, 
and, finally, those satellites which were only included in BR1 la,b due to their post-2004 discovery. Three of 
dwarf initially identifi ed as fossils show in teresting properties (1) a small fraction of the stars in Draco are 
of intermediate age (ICioni & Habingll2005l) . (2) Leo T has ~ 1O 5 M of gas and a yo ung stellar populatio n 
dde Jong et al.ll2008h . and (3) Cetus shows evidence for star formation through z ~ 1 (IMonelli et al]|2010f) . 
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Though we are agnostic about its status, Segue 1 may be an exception. Recent work (ISimon et al. 



Martinez et al.ll2010h suggest that its stellar population has remained well within its tidal radius (thus tides 



are not important) and its stars are unaffected by interactions with the Milky Way. However, other work 



sugge sts that Segue 1 is a highly disrupted star cluster or dwarf (INiederste-Qstholt et all 2009; N orris et al. 
2010). We note that, if Segue 1 is a undisrupted dwarf, the high concentration which has protected Segue 
Fs stars also identifies it as a rare object formed in a high sigma peak at high redshift. The 1 Mpc 3 volume 
of our pre-reionization simulations does not represent a large enough volume to contain a Segue 1. If Segue 
1 is an undisrupted dwarf, than yes, if there are more than one or two additional Segue 1 like objects in the 
Milky Way halo it is a problem for our model that produces larger half-light radii than Segue l's. However, 
if Segue 1 is disrupting than (i) we would not expect to see objects of that type beyond 100 kpc from the 
Milky Way, and (ii) the presence of additional Segue 1 objects would not pose a problem.s 



The ultra-faint dwarfs in the first group, and the classical dSph noted in iRicotti & Gnedinl ((2005) are 
the best candidates for an observed population of primordial fossils, with stellar spheroids not significantly 
modified by tides. However, as noted in Paper I, classical dwarfs with Ly > 10 6 L Q are too bright to be 
hosted in halos with v max < Vfm for Vfut = 20 km s -1 or 30 km s _1 . While they may have formed most 
of their stars before reionization, we exclude them from out comparison to be as conservative as possible. 

Throughout this work, and in Paper I, we compare the observed Milky Way satellites to our luminous 
z = halos with Sy > 10~ L4 L Q pc -2 . We are also able to use our simulations to study the distribution of 
a hereto undetected population of ultra-faints with Sy < 10 -1 ' 4 L Q pc~ 2 and Ly 10 4 L Q . The possible 
existence and undetectability of this population was first noticed in BR09, from the a nalysis of RG05 sim - 
ulations (see also Ricotti 2010 for a review). However, using independent arguments, iBullock et al.1 (|2010i ) 
have also proposed the existence of this population they refer to as "stealth galaxies." 



4. Results 

In this section, we compare the distributions of non-fossils and true-fossils to the galactocentric radial 
distribution of the observed Milky Way satellites. We first compare the galactocentric radial distributions 
of our simulations to observations. We then make detailed comparisons between the observed cumulative 
luminosity function of the Milky Way satellites and the simulated cumulative luminosity functions of our 
non-fossil and true fossil populations. Note, that our simulated cumulative luminosity functions only include 
stellar populations formed before reionization. Therefore, we refer to our simulated cumulative luminosity 
functions as primordial cumulative luminosity functions. Any star formation that may take place in halos 
with v max > Vfm after reionization is not accounted for in our simulated luminosity functions. Thus, 
only the cumulative luminosity function of true fossils can be directly compared to observations, while the 
luminosities of the non-fossils are lower limits. 
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Fig. 1. — Left Galactocentric radial distribution of all simulated satellites for MW.2 and MW.3 from Run D 
(black curves) compared to the radial distribution of all observed Milky Way satellites (red triangles). We 
have included all simulated subhalos and known satellites regardless of their classificati on or whether they 
are detectable. Right Same as the left panel but we have convolved our populations with lWalsh et al.1 (|2009) 
detection limits and only included simulated subhalos which can be detected in the SDSS data. 



4.1. Radial Distribution of Fossils Near Milky Ways 



Figure Q] shows the galactocentric radial distribution of all the simulated and observed Milky Way 
satellites. In the left panel of F igure Q] we compare observations to si mulations without correcting for the 



sensitivity limits of the SDSS (IWalsh et al.ll2009l : iKoposov et al 



2008) or whether a satellite is a fossil. In 



the right panel, we show all the satellites again, now applying the lWalsh et all ([2009) limits to the simulated 
halos around MW.2 and MW.3. Figure [2] shows the galactocentric radial distribution for only the observed 
and simula ted fossils. As in Figure [TJ the right and left panels show the simulated true fossils with and 
without the ' 



Walsh et al 



(2009) corrections. The observational and theoretical fossil definitions are discussed 
in Sections l3l and |2~T1 respectively. Our simulations do not account for tidal stripping of stars, and do not 
reproduce the properties of the inner ultra-faint dwarfs, and we do not include them in Figure [2] 

The left panels of Figures Q] and [2] show that at Ly ~ 1O 5 L the fossils become a significant fraction 
of the satellite population, with fossil dominance increasing as satellite luminosity decreases. This is further 
illustrated in Figure [4] which shows the fraction of subhalos which are fossils, Nj os /N a a, as a function 
of distance from the host for the same luminosity bins as Figure [2l excepting Ly > 1O 6 L . We find that 
for 10 5 L Q < Ly < 10 6 L Q bin, the fraction of fossils is 0.05-0.1, with the fraction decreasing as host 
halo mass increases. For the lower luminosity bins, Nf os /N a u converges to 40 — 50% and 70 — 80% for 
10 4 - 1O 5 L and 10 3 - 1O 4 L bins, respectively. If we include the inner ultra-faints in our galactocentric 
radial distributions, we find a significant overabundance of observed dwarfs within 50 kpc of the Milky Way. 
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Fig. 2. — Left Same as Fig. Q] but the observed satellite distribut i ons o nly include bona fide fossils: the 
classical dSph which were designated fossils in Ricotti & Gnedinl (12005b and the ultra-faints whose stellar 
properties match those of the simulated fossil population. Note that this excludes most of the ultra-faints 
within 50 kpc. In our simulated distributions we use vjuter = 20 km s _1 to define a fossil. We have 
included all simulated fossils, including those which would sit below the SDSS detection limits. Right 
Sam e as the left panel but simulated radial distributions only include the true fossils which would fall within 
the lwalsh et all (|2009T) detection limits. 



The stellar properties of the inner ultra-faint dwarfs do not agree with the simulated stellar properties of the 
fossils. As discussed in Paper I, and §|3]of this work, we argue the majority of these objects may represent 
a population of tidally stripped remnants of once more luminous dwarfs. A possible exception, Segue 1, is 
discussed in § [3] Our simulations are also unable to reliably resolve z = halos within 50 kpc of the Milky 
Way. We therefore have excluded anything with R < 50 kpc from our comparisons. 

Without the ultra-faints with R < 50 kpc, the right panels of Figures Q] and [2] show good agreement 
between the simulated satellite distributions of the true fossils around MW.2 and MW.3 and the observed 
Milky Way galactocentr i c radi al distribution. When we convolve our simulated satellite populations with the 
limits from lwalsh et al. ( 2009 ). we find that the agreement between the distribution of dwarfs around MW.2 
and MW.3 and that observed around the Milky Way agree at all radii and luminosity bins for Ly < W 6 L & 
(see right panel of Figure [2]). We thus argue that, in addition to matching the stellar properties of the ultra- 
faints, our simulated fossils also agree with their galactocentric radial distribution. 

Figure [3] shows the galactocentric radial distribution of the undete cted fossils in our si mulations, after 
excluding detectable fossils according to the detection criterion from IWalsh et all (120091) . We have not 
included the bins with Ly > 10 5 L Q because there are no undetected fossils in this luminosity range within 
1 Mpc of the Milky Way. In addition, all fossils with Ly > W 4 L & are detected within 200 kpc. For 
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the lowest luminosity fossils (Ly < 1O 4 L ) we find ~ 400 — 500 undetected dwarfs within 1 Mpc and 
150 within 200 kpc. We have included a panel for very low luminosity bin (10 2 — 10 3 L Q ) to look at the 
distribution of the dimmest fossils which are invisible beyond a few tens of kpc. While the shape of the 
distribution in the lowest luminosity bins is similar there are approximately two times fewer undetected 
fossils in the 10 2 — 1O 3 L bin. Given that fewer of the fossils in this bin would be detected compared 
to its higher luminosity counterpart, we are seeing the decline of star formation in the minihalos with the 
lowest mass. There are simply fewer 10 2 — 1O 3 L pre-reionization fossils around the Milky Way than their 
10 3 — 1O 4 L counterparts. 

On the other end of the luminosity spectrum, the right panel of Figure [Qshows that while our simulated 
fossils are able to reproduce the ultr a-faint distribution, we see too many massive, bright (Ly > 1O 4 L ) 
satellites at R > R v i r , even after the IWalsh et al.1 (|2009T) corrections are applied. We note that this discrep- 
ancy does not exist in the lowest, 10 3 — 1O 4 L luminosity bin. This is the first evidence of an apparent 
discrepancy between simulations and observations we refer to as the "bright satellite problem." In the next 
sections, we will analyze this discrepancy, try to understand its origin, and whether it can be removed while 
maintaining the agreement of the simulations with observations at smaller radii and lower luminosities. 



4.2. Primordial Cumulative Luminosity Functions 

We next explore the fossil distribution and "bright satellite problem" from another angle via compar- 
isons between simulated cumulative primordial luminosity functions and the observed cumulative luminos- 
ity function at different galactocentric distances from the Milky Way center. Results are equivalent for all 
three simulated Milky Ways and for both versions of our initial conditions. Therefore, for the remainder of 
this section and the next we will be discussing the results for MW.3 in Run D. 

Each cumulative luminosity function in this paper is split into four radial bins to probe different 
regimes. We choose not to include the sample at R < 50 kpc, where the observational sample is the 
most complete, because tidal effects are prevalent and our simulations do not have sufficient resolution to 
determine whether or not pre-reionization halos stripped of their enveloping cloud of tracer particles have 
been tidally disrupted. The first bin we consider shows 50 kpc < R < 100 kpc. The next bin out, the outer 
portion of the Milky Way halo from 100 kpc to 200 kpc, has observations which are fairly complete for 
Ly > 10 4 L , including the brightest ultra-faints. From 200 kpc to 500 kpc all but two of the ultra-faints 
(CVnl and Leo T) would be below the detection limits of the surveys and would not be visible. Roughly, 
this region corresponds to the virial radius (-R200 ~ 200 kpc) to R^q for a Milky Way mass halo. Specific 
to our Local Group, this is the regime where M31 begins to play a significant role in the satellite counts, 
increasing the care required to separate the Milky Way and Local Group dwarfs from those bound to M31. 
The final radial bin, from 500 kpc to 1 Mpc, probes the transition region from the edge of the Milky Way 
halo to the surrounding filament and void. Subhalos at these radii are just beginning to fall into the host 
system, and all the ultra-faints are below detection limits. 

We divide the simulated satellites into fossils and non-fossils: Figure[5]shows the primordial cumulative 
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luminosity functions in the four radial bins for fossil thresholds vtm = 20 km/s (left panel) and Vfm = 
30 km/s (right panel). The observed cumulative luminosity function is shown as magenta lines and includes 
all the classical dwarfs and the ultra-faints, excepting the population at R < 50 kpc. The simulated non- 
fossils are shown as the red solid curve, and for all bins they dominate for Ly > 10 4 — 10 5 L Q . These 
halos may have been able to accrete gas and form stars after reionization, and their primordial cumulative 
luminosities represent a lower limit for their present day luminosity. If we were to allow for additional 
star formation after reionization, the total number of luminous non-fossils would remain constant, but the 
curve would shift to higher luminosities (to the right). The primordial cumulative luminosity function of 
the true fossils (blue dashed curve) has no such caveat. Their luminosities are known since they have 
not undergone post-reionization baryonic evolution aside from the aging of their stellar populations. The 
primordial cumulative luminosity function of the entire simulated population is the solid black curve. Note, 
that for Ly < 10 5 L Q , the total primordial cumulative luminosity function is increasingly dominated by 
fossils. 

Before looking at the primordial cumulative luminosity functions in detail, we insure we are comparing 
equivalent populations. By definition, all observed Milky Way satellites are above current detection limits, 
however, as seen in Paper I, a subset of our simu l ated fo ssils have surface b r ightne sses below the detection 
limit of the SDSS. We use both the Walsh et all d2009l) and iKoposov et all 020081) limits to test the distri- 
bution of detectabl e true fossils again st observations. Figure [6] shows the true fossil luminosity function 
convolved with the IWalsh et al.1 (120091) limits (blue dashed line) and the observed fossil sample (magenta 
line) used in Figured As in the galactocentric radial distributions, we find good agreement between the pri- 
mordial luminosity function of the true fossils and the observed fossils for 50 — 100 kpc and 100 — 200 kpc. 
We do not make comparisons at R > 200 kpc because of the inability of current s urveys to detect fossils 
at these larger distances. In Figure|71 we next use the surface brightness limits from IKoposov et all ((2008) 
to remove any simulated fossil satellite not detectable by current surveys. Using the Koposov et al straight 
surface brightness cut all but eliminates the fossil population for vjut = 20 km s _1 . This is a much stronger 
effect on the detectability of our true fossils than that seen for the Walsh et al. luminosity and distance cuts. 

In all distance bins, there is an overabundance of the bright satellites at luminosities t ypical of the 
classi cal dwarfs (Ly > 10 5 L Q ). These should be easily detectable by the SDSS according to IWalsh et al. 
(|2009r ) (assuming the undetected dwarfs have the same distribution of half light radii as the ultra-faints). 
In Figure [5J the detectable dwarfs are to the right of the dashed line. During our discussion of the missing 
bright satellites, we use cumulative luminosity functions which have not been corrected for the SDSS limits. 
We now look at each distance bin individually. 



4.2.1. Inner Ultra-Faint Dwarfs at R < 50 kpc 

In Paper I, we argue that, while the inner ultra-faints have likely lost significant fractions of their stellar 
populations to tidal striping, they were not necessarily dlrr at the start of their encounters with the Milky 
Way. Instead, they may have been more massive primordial fossils. We base this conjecture on Figure 12 in 
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Paper I which shows that the inner ultra-faints have metallicities, [Fe/H], that are similar to fossils that are 
slightly more luminous. But are there enough massive fossils to account for the inner ultra-faints? Figured] 
shows the mass function (bottom) and luminosity function (top) of the pre-reionization halos which are not 
part of a bound halo at z = and are between 20 kpc and 50 kpc from MW.3. The dotted horizontal lines 
show the approximate number of stripped fossils required to reproduce the inner ultra-faints. We see that 
to produce the ~ 30 inner ultra-faints around the Milky Way, we would only need to consider the largest 
primordial fossils with masses at reionization M > 10 s M and initial luminosities Ly > 10 6 L . 

4.2.2. 50 kpc<R< 100 kpc 

A strong piece of evidence for the primordial model would be the total number of observed satellites 
in one or more radial bins being greater than the number of non-fossils. When we look at R < 100 kpc 
without the R < 50 kpc cut, we see such an overabundance of observed dwarfs. However, when we do not 
include the dwarfs within 50 kpc of the galactic center the case is no longer clear cut. If the satellite count 
from 50 — 100 kpc increases to greater than 25, there is a case for fossils even using the most conservative 
Vfut = 20 km s _1 . For Vfm = 30 km s _1 the number of non-fossils available from 50 — 100 kpc drops to 
~ 18. 

For luminosities at which the observational sample is complete, to the right of the dashed lines, we see 
too many bright (Ly > 10 4 L Q ) objects, even in this inner most radial bin. In addition, as our simulations 
do not account for post-reionization star formation, it is likely that the overabundance of bright objects is 
worse than shown in Figure [5] Unless all of the non-fossils have accreted no gas and formed no additional 
stars after reionization, the simulated curve must lie below the observations. This allows these star forming 
halos to increase in luminosity, shifting the luminosity function of the non-fossils to the right. 

4.2.3. 100 kpc < R < 200 kpc 

In this bin, we probe the outer reaches of the Milky Way's virial halo and there are a few notable 
characteristics of the luminosity functions. First, with the addition of the observed dwarfs in this bin, the 
total number of known satellites around the Milky Way increases to ~ 65, ~ 45 not including the inner 
ultra-faints. The sample at these larger radii is only complete for Ly > 10 4 L Q . For 100 — 200 kpc, with 
a Vfut = 20 km s _1 and the less conservative Vfm = 30 km s _1 there are ~ 60 and ~ 40 simulated non- 
fossils respectively. Second, the presence of undetected dwarfs is corroborated by the shape of the observed 
luminosity function around 10 4 L Q . Not only is it rising steeply to the detection limit at R = 100 kpc, but 
its shape is similar to the simulated primordial luminosity function for the true fossils. 

In the outer virial halo, we once again overproduce the number of bright satellites. At these radii the 
discrepancy between theory and observation is more severe than for 50 kpc< R < 100 kpc since for these 
radii there is only one observed satellite with Ly > 3 x 10 4 L Q . 
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4.2.4. R > 200 kpc 



For R > 200 kpc, we can only make observational comparisons for Ly > 10 5 L Q . Beyond the virial 
radius the discrepancy between the observed number of satellites and our simulations is up to ~ 1.5 orders 
of magnitude, compared with factors of ~ 2 and ~ 10 for the 50 kpc < R < 100 kpc and 100 kpc< R < 
200 kpc bins, respectively. 



Our simulations of the fossils of the first galaxies are consistent with the observed Milky Way satellite 
galactocentric radial distributions (Section |4~TT) . primordial cumulative luminosity functions (Figure 0, and 
internal stellar properties (see Paper I). More intriguing, is that at first glance the model appears to fail in the 
outer parts of the Milky Way by over-producing the number of bright non-fossil satellites. In this section, 
we explore possible solutions to the bright satellite problem and whether the proposed solutions maintain 
the agreement between the observed ultra-faint dwarf population and our simulated true fossils. First, we 
will explore whether the model overrestimates the star formation efficiency in pre-reionization dwarfs, and, 
then, whether we have too many luminous galaxies forming before reionization. It would be of great interest 
if we could use current observations to constrain galaxy formation models before reionization. We conclude 
the analysis with a proposal in which the bright satellites do exist in the outer parts of the Milky Way halo 
but may still be elusive to detection due to their extremely low surface brightnesses. The discovery of 
these "ghost halos" around known isolated dwarfs is a test of the primordial model and would allow us to 
determine the efficiency of star formation in the first galaxies. 



First, we explore whether our pre-reionization simulations overestimate the luminosity of galaxies in- 
dependently of dark matter halo masses. This could be due to using an incorrect IMF in the pre-reionization 
simulations. We use a mass-to-light ratio for the aged stellar population, 



where M^ ei and is the mass of the stellar population at reionization, and M* y is the mass of the stellar 
population at z = 0. The ratio of Ml ei /Mt° day depends on the primordial IMF. The ratio of M^ ei /Mi° day 
depends on the primordial IMF, and M l ° day jL 1 — 2 as for the oldest globular clusters in the Milky Way. 
A more top heavy IMF for the primordial stars will result in greater stellar mass loss after reionization and 
fewer low mass stars which can survive to the modern epoch. Conversely, the ratio drops as the primordial 
IMF produces fewer high mass stars. 

To approximate this effect, in Figure |9j we plot the cumulative luminosity functions as in Figure [51 but 



5. Where are the Bright Satellites? 



5.1. Increasing Mass to Light Ratios 




(1) 
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increase the stellar mass to light ratio by a factor of 10 in our pre-reionization dwarfs, to 50Mq/L & . The 
figure shows that increasing the mass to light ratio to 50M & / L & does not decrease the number of luminous 
satellites enough to match observations. In the 50 kpc < R < 100 kpc bin, we can match observations. 
However, since the primordial luminosities of the non-fossils are only lower limits, the agreement disappears 
if the population formed any stars after reionization. If they did, we are still over-producing luminous 
satellites. We need to use a mass to light of 500M & /L & to not over-produce the number of non-fossil 
satellites in any radial bin. However, this high mass to light ratio makes the fossils virtually dark, with 
L v < 10 2 L . 

A blanket suppression of star formation in all halos does not solve the bright satellite problem unless 
we suppress all star formation in most halos before reionizatioin. We next explore suppression mechanisms 
which are dependent on the environment or properties of the halos. 



5.2. Suppression of Pre-Reionization Dwarf Formation in Voids 

The formatio n of H2 in the early universe is catalyzed by ionizing UV radiation emitted by nearby star 



forming regions (IRicotti et al.H2002bl) . In the voids, two factors work against H2 formation. The delay of 
structure formation in the voids relative to higher density regions will prevent the minihalos from collapsing 
until lower redshifts when the H2 dissociating background is stronger. In addition, the importance of positive 
feedback is reduced due to the l arger m ean distances between minihalos in the voids and sources of ionizing 



radiation (Ric otti et al.ll2002alJbLl2008r) . The combination of these factors may result in a reduced abundance 
of H2 relative to the regions around a Milky Way. This may produce a star formation efficiency before 
reionization that depends on the environment. We approximate the most extreme case of H2 suppression in 
the voids by suppressing all star formation in halos in regions with 5 0.4. The extreme suppression of the 
star formation in the voids we use treats all halos but those in the overdense regions (z e ff = 8.3) as dark. 

Since the bright satellite problem is most prominent in the outer regions of the Milky Way halo, the 
lack of star formation in low density regions may decrease the number of bright halos beyond the virial 
radius while leaving the satellite luminosity functions unchanged at smaller radii and lower luminosities. 
However, Figure [H] shows that, even in the most extreme case, suppressing H2 formation in the voids does 
not decrease the number of bright satellites enough in any radial bin to bring our simulations into agreement 
with observations. In this scenario, there is no appreciable reduction of the number of bright satellites for 
R < 200 kpc. There is a decrease in the luminosity function at larger radii, but it is neither strong or 
focused enough on the high luminosity subhalos to solve the overabundance of bright satellites in the outer 
parts of the Milky Way. Enough of the region within 1 Mpc of our Milky Way is at the highest density in 
our simulation (z e ff = 8.3) that the complete suppression of star formation in even moderately less dense 
regions does not sufficiently change the luminosity functions or radial distributions. 



With the inability of H2 suppression in the voids and overall suppression of star formation to account 
for the missing bright satellites, we shift our focus to properties of the halos which vary with halo mass. 
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5.3. Lowering the Star Forming Efficiency 

In the pre-reionization simulations, the sub-grid recipe for star formation depends on a free parameter, 
e*, controlling the efficiency of conversion of gas into stars per unit dynamical time. One of the main 
results of the pre-reionization simulations is that the global star formation rate and /*(M) = M*/Mb ar 
(the fraction of total mass converted into stars, where M& ar = M/7.5) is nearly independent of e* in 
small mass dwarfs due to the self-regulation mechanisms of star formation. However, in halos with masses 
M 5 — 10 x 10 7 Mq, f*(M) is typically proportional to e* since the higher mass minihalos are less 
sensitive to self -r egulation feedback. We used e* = 5% in our fiducial runs, that may be too large (e.g., 



Trenti et al.ll201fj) . The pre-reionization simulations may have overestimated the luminosity of primordial 
dwarfs with M ^ 5 x 10 7 M Q , for which the /* vs M relationship is tighter. We explore the effect of 
reducing e* by introducing a maximum stellar fraction, f* jCr u. Roughly, /* iCr jt corresponds to the mass 
threshold where feedback effects no longer dominate and where the value of e* becomes important. If we 
reduce f* jCr u, we will decrease the luminosities of our most luminous halos. Roughly, for halo masses 
M ~ 3 x 10 7 M (virial mass at formation) our simulations have f* (M ) ~ 1%. This is in agreement with 



observed values for dwarfs with v c ~ 10 km s 1 ( McGaugh et al. 201dh 



Figure[[2]shows the luminosity functions of our simulations with all halos with f* >cr it = 1% (left panel) 
and f* >cr it = 0.1% (right panel). The figure shows that lowering the star formation efficiency preferentially 
for the higher mass halos is effective in decreasing the number of non-fossil subhalos with Ly > 10 5 L Q . 
Adopting /* )Cr it = 1% decreases the number of luminous halos enough to bring the luminosity functions 
in agreement with observations, while preserving the agreement for the fossil population. However, in 
the radial bins 50 kpc < R < 500 kpc, there still are too many subhalos with Ly > 10 5 L , though 
the discrepancy has dropped significantly. Coupling a lower e* with a higher mass to light ratio or H2 
suppression in the voids does not correct the remaining bright satellite overabundance. However, when 
we set /*,crit = 0.1%, the cumulative primordial luminosity function of our simulated dwarfs becomes 
consistent with observations at all radii, but requires a deduction of f* in h alos with mass M ~ 7 x 10 6 M 



whose /* is self-regulated and thus independent of e* (|Ricotti et al.ll2002bl ). 



Any solution for the overabundance of bright satellites must preserve not only the existence of the 
true fossil population, but also its distribution and properties. We next look at the other dimensions of 
the agreement between the true fossil populations with an f cr i t = 1% and 0.1% and the ultra-faints and 
classical dSph. Figure [13] shows the radial distribution of the true fossils around MW.2 and MW.3 from 
Run D the observed Milky Way population for an f cr i t = 1%. While for 1O 3 L < Ly < 10 4 L and 
10 4 L Q < Ly < 10 5 Lq the fossil population with f cr u = 1% reproduces the observed radial distribution, 
we no longer have any true fossils with Ly > 1O 5 L . If the star formation efficiency of our pre-reionization 
halos is lowered enough to bring the number of luminous satellites in line with observations, the fossil 
luminosity threshold discussed in Paper I is dropped to Ly < 10 5 Lq. For f cr u = 0.1% the threshold drops 
further to 1O 4 L . 



The loss of the multi-dimensional agreement between our true fossils and the ultra-faints shows that 
lowering e* enough to account for the missing luminous satellites is not a viable solution for the bright 
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satellite problem if the ultra-faint dwarfs are fossils of the first galaxies. In this interpretation, we need 
a different mechanism which will either preferentially suppress star formation in the most luminous pre- 
reionization halos to a greater degree, or cause their lower redshift counterparts to lose the majority of their 
primordial stellar population after reionization. If neither of these solutions work, we must ask ourselves if 
the halos CDM predicts are there at all. 



5.4. The Ghost Halos 

As discussed in Paper I, our N-body method does not allow us to determine the dynamics of the stars 
in halos that undergo mergers, or the degree to which those stars are tidally stripped. However, we have 
used analytic relationships to estimate the importance of dynamical heating of the stars when z = fossils 
(about 20%) are produced by mergers of more than one pre-reionization dwarf (Paper I) We refer to such 
interactions as galaxy mergers. Here, we focus on those dynamical processes in non-fossils which result 
in the dispersion of the primordial stellar populations of the brightest satellites, the net effect of which is 
to either make the non-fossil populations invisible to current surveys by reducing their surface brightnesses 
below the SDSS detection limits or preferentially stripping them during interactions with more massive 
halos. The former mechanism would be relevant to non-fossils at R ^ 500 kpc where tidal forces are 
negligible. 

The number of pre-reionization halos in a z = dwarf increases with mass. In this section, we explore 
the role of mergers to rend invisible, or strip, the primordial populations of stars in the the more massive 
dwarfs (non-fossils). Unlike in the previous sections, here we differentiate between the non-fossils and 
polluted fossils in our simulations. We remind the reader, that though both populations have v max which 
were large enough for them to accrete gas from the IGM in the past, only the non-fossils are at or above that 
threshold at z = 0. 

When a system undergoes a galaxy merger, kinetic energy from the collision is imparted to the stars. 
Immediately after the collision, the new system will be in its most diffuse state. We define a galaxy merger 
as the interaction of two or more pre-reionization halos, both containing a primordial stellar population. 
Although there is significant scatter in the luminosities of minihalos of the same mass, in general, the 
luminous minihalos are more massive than those which are dark. An interaction between two luminous 
minihalos is therefore more significant. We use N\ um , the number of luminous pre-reionization halos within 
at z = halo, as a proxy for the number of galaxy mergers. If there are multiple galaxy mergers in a short 
amount of time, these stars will be susceptible to stripping. In addition, recent work on increasing the extent 
of the stellar population in bright ellipticals from z = 2 to z = suggests that many minor interactions 
over several Gyrs can increase the size of the galaxy by a factor of 2-5 without significantly increasing mass 
( Naab et aDboOSO . Roughly, the larger the number of significant interactions, the greater the spatial extent of 



the pre-reionization population and more likely the halo will have lost a significant fraction of its primordial 
population to dynamical heating. 

For isolated halos, with large Ni um , the radius of the primordial population increases, possibly until 
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it fills the spatial extent of the dark matter halo. Such an extended system would have extremely low 
surface brightness and would be susceptible to tidal stripping. We next look at which of our three subhalo 
populations has a significant number of members with Ni um > 3. The fractions of the true fossil, and 
non-fossils populations with a given N\ um are shown in the left and right panels of Figure [14] Neither of 
the fossil populations has a significant fraction of subhalos with Ni um > 3 with the fractions at ~ 1% and 
~ 10% respectively. We therefore assume, that, while a few of our true and polluted fossils may have had 
their primordial populations diffused by mergers, the vast majority remain dynamically cold. 

The non-fossils show the opposite trend. The right panel of Figure [14] shows that < 10% of the non- 
fossils have Ni um < 3 and the distribution peaks at Ni um ~ 5. A population of non-fossils would be 
much more likely to have a primordial population dispersed by multiple major interactions than their fossil 
counterparts. 

The non-fossil populations (see Figure [5]) could have either lost their stellar populations, or had their 
primordial populations increase in size to the point where they are undetectable by the SDSS. Near the Milky 
Way, we assume they lost their entire primordial stellar population to tidal interactions. After falling into 
the Milky Way halo, the non-fossils were unable to form a significant younger stellar population. 

At larger radii, the non-fossils are less likely to have their primordial populations stripped. However, 
as the stars expand to fill the spatial extent of the dark matter halo, the non-fossils end up with a primordial 
population with Sy ~ Ly x R p 2 , where Ly is the luminosity of the primordial population and R p is 
the radius of the primordial population. Figure [15] shows the fraction of the non-fossils with a primordial 
population with surface brightness, Sy, for R p = R m ax and R p = 0.25 x R max , where R m ax is the radius 
of the maximum circular velocity. We find that for R p = R ma x about 1% of non-fossils would have an 
extended primordial halo above the SDSS detection limits (to the right of the dashed line), and, when R p is 
decreased to 0.25 x Rmax, the detectable fraction only rises to 20%. If these non-fossils formed few or no 
stars after reionization the majority would be undetectable by SDSS. The non-fossils which did form stars 
after reionization today coulds be dlrr or one of the few isolated dSphs or dSphs/Irrs: e.g., Cetus, Tucana, 
Antlia. In our scenario, they would be surrounded by "ghost halos" of primordial stars ~ 12 Gyr old with 
[Fe/H] £j — 2. But does the dispersal of the primordial population into "ghost halos" solve the "bright 
satellite problem"? 

We quantitatively approximate this for our simulations by using a circular velocity cut. We look at the 
primordial luminosity functions as if all the ghost halos are either stripped or below SDSS detection limits. 
Practically, we set the luminosities to zero for all the non-fossils. We look to see if this cut solves the bright 
satellite problem while preserving the fossils better than lowering the star formation efficiency. We find it to 
be a good solution to the bright satellite problem. 

Since setting the non-fossils luminosities to zero is able to decrease the number of luminous satellites, 
we look at the luminosity function it produces in more detail. All the curves in Figure [16] are the same as in 
Figure [5] excepting the red curve for the non-fossils that now represents only the polluted fossils. 

We now look at each distance bin to see what turning off the non-fossil population has done to our 
various arguments. For 50 kpc < R < 100 kpc, the necessity of a primordial dwarf population is even 
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stronger when we only consider the fossil and polluted fossil populations. There are only ~ 7 polluted 
fossils within 100 kpc, less than one-fifth of what is required to account for the ~ 50 observed satellites. For 
Ly > 1O 4 L , the luminosity function now sits below observations. This gives the remaining star forming 
halos room to form additional stars without overproducing subhalos with Ly > 1O 4 L . The total number 
of subhalos within 100 kpc decreasing to 35 is consistent with observations. First, our MW.3 is on the low 
end of the mass range for the Milky Way for both observational estimates and simulations. Second, as seen 
in Figured! there are more than enough stripped down fossils which formed in halos with M > 1O 7 M and 
Ly > 10 5 L Q to fill in the deficit. Since these objects would lack a cloud of tracer particles, they are marked 
as unbound by the halo finder AHF and would be included in the luminosity and mass of the host halo and 
not in any luminosity function of the satellites. 

The complete invisibility of the non-fossils is not quite as successful for 100 kpc< R < 200 kpc as we 
are still slightly overproducing the number of Ly > 10 5 L satellites compared to observations. However, 
we are better able to reproduce the sudden steepening in the observed luminosity function in Figure [T6l than 
with any of our other suppression mechanisms (Figures [lOl \TT\ and especially IT2T> . This feature may be 
unique to the Milky Way so we are not unduly concerned with matching it. In addition, if the non-fossils 
are dark, our argument for the existence of primordial fossils becomes straightforward. There are only 
~ 30 polluted fossils in this distance bin, only 75% of the ~ 40 observed galaxies, with any dwarf with 
Ly < 10 4 L difficult, if not impossible, to detect with current surveys. 

Beyond the MW.3 virial radius (R ~ 200 kpc), turning the non-fossils dark easily places the primordial 
luminosity function into agreement with observations. This allows for the formation of post-reionization 
populations of stars in the polluted fossils and non-fossils. We remind the reader that the z = halos at 
these radii would likely be on first approach to the Milky Way system and more likely to accrete and retain 
gas at later times. The diffuse primordial population in these distant non-fossils is an observational test of 
star formation in pre-reionization dwarfs and the existence of pre-reionization fossils. 



6. Three Observational Tests and Model Predictions 

Although the existence of pre-reionization fossils seems likely, observations do not unequivocally 
demonstrate their existence due to the large uncertainties in estimating the number of yet undiscovered 



ultra-faint dwarfs jTollerud et al. 2008 ). In this section, we summarize three observational tests for the ex- 
istence of fossils of the first galaxies that we propose based on the results in Paper I and the present work. 
The first test of our model is especially interesting as it can be performed using HST observations and does 
not require waiting for future all sky surveys deeper than SDSS, like PanStar or LSST, to be online. 

1. "Ghost halos" around dwarfs on the outskirts of the Local Group 

The primordial stellar populations in minihalos that formed before reionization should produce diffuse 
"ghost halos" of primordial stars around isolated dwarfs. We have shown that the total luminosity of 
the "ghost halos" is comparable to the one of classical dwarfs, but the surface brightness of the stars 
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is well below the SDSS detection limits. Contrary to the difficulties of rinding ultra-faints, we know 
where these diffuse stars are and we can plan deep observations to detect them. The diffuse primordial 
stellar populations around non-fossils should not be tidally stripped in dwarfs with galactocentric 
distance > 1 Mpc from the Milky Way. We do not know within which distance tidal stripping would 
become important, but due to their large half light radii, they certainly are the stellar population 
that would be stripped first. "Ghost halos" can be best detected by resolving their individual main 
sequence stars around isolated dlrrs or dSphs before using spectra to determine their metallicities and 
dynamics. Unlike younger stars dispersed from the central galaxy, the primordial ghost halo would 
have a [Fe/H] < —2.5, and we are currently running simulations to determine their dynamics. 

Recent HST observations of M3 1 have resolved the main sequence using ACS (brown et al. 2006 , 



2008, 2009|). A low-luminosity dwarf at ~ 800 kpc on the other side of the Milky Way from M31 



would be an excellent candidate for the ghost halo search, and detection of its main-sequence pri- 
mordial stars would be within the reach of HST. At 800 kpc, the field of view of WFC3 (162") is 
~ 630 pc, at 1 Mpc, ~ 785 pc, and at 2 Mpc, ~ 1.6 kpc. The ghost halos are > 1 kpc in radius, 
often up to a few tens of kpc, therefore, WFC3 would not be able to image the entire dwarf with its 
ghost halo. However, aimed at the outskirts of a likely ghost halo host, it could look for signs of a 
primordial halo in the color magnitude diagram and radial surface brightness distribution. We would 
be able to resolve the individual stars in the ghost halos at ~ 1 Mpc. The ghost halos have surface 
densities of stars of 0.001 and 1 star pc~ 2 depending on the extent of the ghost halos and the slope 
of the IMF at low masses. Assuming a stellar density of 1 star pc~ 2 , the angular distance between 
each star is ~ 0.26", larger than the WFC3 resolution of 0.04" per pixel. Since the stellar population 
can be resolved, determining the details of a ghost halo population is a matter of taking deep enough 
exposure to detect the main sequence stars. Red giant branch stars, while brighter and easier to detect, 
have a density three orders of magnitude lower than the main sequence. 

Deep observations work well for detecting the ghost halos when we already know where they are. 
These primordial populations surround dwarf galaxies that have undergone significant star formation 
since reionization and may have detectable gas and active star formation today. However, as discussed 
in § 15 -41 in order to reproduce the observed satellite distribution, only a fraction of the ghost halos can 
have formed a significant younger stellar populations. We remind the reader that our definition of 
a fossil versus a non-fossil in the simulations assumes a constant Vfut er = 20 km s _1 . Assuming 
a larger filtering velocity produces a smaller number of non-fossils, some of the halos we define as 
non-fossils, and containing ghost halos, may have been unable to accrete gas and form stars after 
reionization due to additional heating of the IGM. This population of "dead" ghost halos would have 
only an extremely diffuse, primordial population. Without H I or more concentrated, younger stars, 
the best chance for detecting these ghost halos would be large scale surveys. Figure [15] shows that 
the majority of the ghost halos are beyond the reach of the SDSS, but what about upcoming, deeper 
surveys such as PanSTARRS? 

Diffuse stellar systems like the ultra-faint dwarfs and the ghost halos are found in surveys by looking 
for overdensities of stars relative to the background. In many cases, by detecting stars at the tip of the 
red giant branch (RGB). In low luminosity systems the detection of the RGB depends on two factors, 
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the distance to the halo and the population of the RGB. Low luminosity systems, like the ultra-faint 
dwarfs, can have as few as a thousand stars, and therefore a sparsely populated, and difficult to detect, 
RGB. For example, with SDSS (magnitude limit r = 22.5) Hercules (1.1 x 1O 4 L ) could only be 
detected to 300 kpc, while the mor e luminous CVn I (2.3 x 1O 5 L ) would be seen at a Mpc from 
the Milky Way ( Ko posov et al. l l2Q08h . PanSTARRS (magnitude limit r=24) will reach 1.5 magnitudes 



deeper that SDSS, detecting the same RGB twice as far. However, the primordial populations in 
the ghost halos are extremely diffuse, and it is unclear that the overdensity of their RGB stars would 
be high enough to be detected against foreground M dwarfs and distant galaxies. In short, while 
PanSTARRS is expected to detect new ultra-faint dwarfs, it may not the best tool for finding ghost 
halos. 

If some of the discovered ultra-faint dwarfs are fossils, then "ghost halos" should exist. Vice versa, 
the detection of "ghost halos," regardless of the method, can be used to constrain the star formation 
rates before reionization and would imply the existence of fossils, although these fossils may not have 
yet been discovered due to their low surface brightnesses. 

2. Population of yet undetected ultra-faint dwarfs 

In BR09 and Paper I we have discussed in detail the existence and the properties of a yet undetected 
population of fossils with Ly < 10 4 L and surface brightness Sy < 10 -1 ' 4 L pc~ 2 . This popu- 
lation of ultra-faints should be accessible to future all sky survery such as PanStars and LSST If the 
undetected dwarfs are fossils, they should have stellar velocity dispersions equivalent to those of the 
ultra-faints with corresponding mass to light ratios of 1O 4 M /L or higher. The typical [Fe/H] of 
the "stealth" fossils should be < —2.5. We direct the reader to Paper I for a detailed justification of 
these predictions. 

3. Dark and ultra-faint gas rich dwarfs in the voids 



According to the model proposed in lRicot ti (2009), a subset of minihalos in the voids may have been 



able to condense gas from the IGM after Helium II reionization (at z ~ 3). However, they would not 
form stars unless their gas reached a sufficient density. These minihalos may or may not have formed 
stars before reionization, and any stellar populations they did have would be below the detection limits 
of both current and future surveys. The HI in these objects could be detected by blind 21 cm surveys. 
Recently, ALFALFA and GALFA surveys have reported the discovery of several small and compact 
clouds of neutral hydrogen, some of which may represent a population of pre-reionization minihalos. 
Some of these clouds could be "dark galaxies" and represent the smallest detectable halos around the 
Milky Way, others may be ultra-faint dwarfs in the voids. The location of unassociated HI detections 
could then guide optical surveys to these primordial fossils in a focused deep search for their ancient 
populations. 
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7. Summary and Discussion 

Through this work and Paper I, we have used results of simulations to study the origin of observed 
Milky Way and M31 satellites and understand whether they are compatible with models of star formation 
before reionization. In the primordial model, a subset of the Milky Way satellites formed with their cur- 
rent properties with minimal modifications by tidal stripping. These low luminosity satellites formed the 
majority (> 70%) of their stars in minihalos before reionization. Our simulated true fossils produce an 
excellent agreement in properties and distribution with the observed lowest luminosity Milky Way satellites. 
In BR09 and Paper I, we showed that a subset of the ultra-faints, all with Ly < 10 5 L , have half-light radii, 
surface brightnesses, mass-to-light ratios, velocity dispersions, metallicities and metallicity dispersions con- 
sistent with the expected stellar properties of the true fossils. In the present paper, we have compared the 
galactocentric radial distributions and primordial cumulative luminosity functions of simulated fossils to 
observations of Milky Way satellites. When we compare the observed and simulated distributions, we find 
them to be in agreement with each other for Ly < 10 5 L Q . In addition, a large population of primordial 
fossils have surface brightness below the detection limits of current surveys (e.g., SDSS). The following list 
summarizes the main results of this second paper of the series. 

• We are able to reproduce the distribution of the ultra-faints with our simulated primordial fossils. We 
find no missing satellites at the lower end of the mass and luminosity functions, but our model predicts 
~ 150 additional Milky Way satellites detectable by upcoming surveys (PanStars, LSST). In addition 
to the undiscovered "missing" satellites at distances too great 50-100 kpc) to be detected by SDSS 
or in areas on the sky not in the SDSS footprint, there is also a new population of ultra-faints with 
mega-faint surface brightness hidden below the SDSS detection limits at all distances from the Milky 
Way. 

• At all radii, we find an overabundance of simulated bright (Ly > 1O 4 L0) satellites with respect to 
observations, which, even with only their primordial luminosities, would be easily detected by current 
surveys if we assume they have the same half-light radii as the other known ultra-faint dwarfs. Given 
the agreement between the stellar properties and distributions of the ultra-faints and those of our fossil 
dwarfs, we cannot account for the excess bright satellites by imposing a blanket suppression of star 
formation below a given mass. 

• Lower H2 formation rates and subsequent lower minihalo star formation rates in the voids are not able 
to bring the number of bright satellites into agreement with observations. 

• Effectively lowering the star formation efficiency can fix the bright satellite problem if we assume 
only pre-reionization halos with M < 7 x 10 6 M Q had SFR dominated by local, stochastic feed- 
back. However, not only is this contrary to current understanding of star formation in minihalos, 
but the fossil population this "solution" produces cannot reproduce the distribution of the ultra-faint 
population. 

• We bring the number of bright satellites into agreement with observations, while leaving the fossil 
population untouched, by assuming the primordial stellar populations of our non-fossils (with max- 
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imum circular velocities, v max (z = 0) > vju t ) become extremely diffuse via kinetic energy from 
galaxy mergers. The existence of "ghost halos" of primordial stars is a new powerful observational 
prediction of our model that can be straight forwardly tested using HST observations of isolated 
dwarfs around the Local Group. 



One of the key predictions of the primordial model is a total number of satellites for the Milky Way 
between 200 — 300, only a maximum of 100 of which are non-fossils (here we assumed Vfu t = 20 km s _1 ). 
The number of Milky Way satellites which are not fossils provides an important test for star formation in 
minihalos at high redshift. If, after PanSTARRS and LSST are online, the number of ultra-faint Milky Way 
satellites remains < 100, we have a strong constraint on star formation in pre-reionization dwarfs. Either no 
pre-reionization fossils survived near the Milky Way, or almost none of the halos with masses at formation 
M < 1O 8 M formed stars. However, if the satellite count rises to > 100, some of the dimmest Milky 
Way dwarfs must be fossils of reionization. Using details of the stellar populations and their distributions, 
observations of these fossils can constrain models of star formation at high redshift. 

A caveat to this picture is that the number of non-fossils is highly sensitive to the choice of the filtering 
velocity. When we raise the filtering velocity to 30 km s _1 , the number of non-fossils drops by a third to 
60±8 from the 90±10 for Vfut = 20 km s _1 . The choice of 20 km s _1 assumes a constant IGM density with 
Tjqm = 10 4 K throughout a minihalo's evolution. In reality, the situation is not so simple. The gas near 
10 12 M Q halos and in the filaments between may be heated to ~ 10 5 — 10 6 K by AGN feedback. The higher 
temperatures of this local intergalactic medium may correspond to vtm £ 40 km s _1 . In addition to the 
higher filtering velocity, the higher density near a Milky Way mass halo reduces the effective potential depth 
of the subhalos, increasing the mass threshold for post-reionization gas accretion still further. Simulations 
to determine the temperature and density of the IGM near a Milky Way from reionization to the modern 
epoch are needed to determine the (x, z), and whether these factors can explain the existence of Milky 
Way and M31 dwarfs with the observed properties of fossils, but luminosities above the 10 6 L threshold. 

The observed distributions we compare to depend on how we correct for the incomplete sky coverage 
of the SDSS. In this work, we have assumed an isotropic satellite distribution at R > 50 kpc. Under this 
assumption, the SDSS completeness correction for the ultra-faints is 3.54. We briefly check if the agreement 
between the observed a nd simulated distributions is dependent on the isotropic assumption. Recent work 



(IMetz et al.1120 07. 2009) has suggested that rather than being isotropic, the Milky Way satellites are oriented 
in a plane approximately perpendicular to the disk. We approximate this non-homogeneous satellite distri- 
bution by correcting for the SDSS sky coverage by a factor of 2.0 instead of 3.54. The number of classical 
fossils remains the same. The different correction does not change the consistency of our simulated galacto- 
centric distribution with observations, though the lower correction factor suggests a higher Milky Way mass. 
It also does not change the bright satellite problem, in fact, the lower observational correction factor makes 
the overabundance of simulated Ly > 10 4 L Q dwarfs worse by about a factor of two. 

We have suggested two solutions which correct for the overabundance of bright satellites while pre- 
serving at least a fraction of the primordial fossil population. Each presents a different picture when we 
consider it in the context of the voids. The first, and less effective, solution calls for a low star formation 
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efficiency. In this picture, the 10 10 M halos visible in current surveys will have their star formation damp- 
ened, however as we move to 10 7 M we enter the regime where stochastic feedback effects dominate over 
the choice of e*. Thus, the voids would appear relatively empty, but only because we cannot yet detect the 
less than 10 5 L Q fossil populations which formed in the 10' M Q halos before reionization. 

The dispersal of the non-fossils' primordial populations into ghost halos is a more effective solution 
to the "bright satellite problem" within 1 Mpc of the Milky Way, but leaves a conundrum in the voids. 
Regardless of whether the primordial population would be detectable, how do we keep the post-reionization 
star formation in these non-fossils low enough to prevent this later star formation from producing more 
My > — 16 galaxies than are currently observed? 

Any post-reionization star formation in the non-fossils results in a young population which would be 
(i) brighter and bluer, (ii) more concentrated, since enriched gas will cool faster and sink deeper into the 
gravitational potential, and (iii) possibly accompanied by an H Ireservoir. Any of these properties would 
make the post-reionization population easier to see, and the the non-fossil harder to hide. To suppress 
the post-reionization baryonic evolution in the non-fossils we examine our naive assumption that they all 
undergo significant baryonic evolution after reionization. 

The easiest way to suppress star formation in the lower mass non-fossils is to raise the filtering velocity. 
As has already been discussed in Chapter 4, any non-fossil embedded in the WHIM (T ~ 10 5 K ) would 



have a Vfut er ~ 40 km s _1 . However, the WHIM does not exist until z < 1 dSmith et alJl2010h . leaving 
~ 6 Gyr after reionization when the non-fossils could have accreted gas and formed stars. Active galactic 
nuclei (AGN) eject enormous amounts of energy into their environs, heating the gas and raising the filtering 
velocity, but how far from the host galaxy the AGN is effective at suppressing star formation in dwarfs ,and 
for how long is unclear. In addition, we see a bright satellite problem around every large Local Volume 
galaxy. Is it reasonable to assume that, at some point in its evolution, every L* galaxy hosted an AGN? A 
final possibility is that reionization was extremely efficient at quenching star formation in 20 — 40 km s" 1 
halos and the non-fossils were never able to build up enough gas from the post-reionization IGM to form 
additional stars. 

In summary, while the bright satellite problem can be "solved" for the primordial population alone, we 
still need to account for the post-reionization evolution of the non-fossils. In order to maintain the agreement 
with observations, only ~ 10% of the non-fossils can form significant stellar populations after reionization. 
Determining how and if the other ~ 90% can be suppressed will tell us how much of a problem the bright 
satellite problem is. 



The simulations presented in this paper were carried out using computing clusters administered by 
the Center for Theory and Computation of the Department of Astronomy at the University of Maryland 
("yorp"), and the Office of Information Technology at the University of Maryland ("hpcc"). This research 
was supported by NASA grants NNX07AH10G and NNX10AH10G. The authors thank the anonymous 
referee for constructive comments and feedback. Thanks from MSB and MR to Stacy McGaugh, Derek 
Richardson and Rosie Wyse for helpful conversations and comments. MSB would like to thank Susan 
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R < 50 kpc 


R > 50 kpc 


Inconsistent 


Segue 1 
Segue 2 
Willman 1 


Pisces II 
Leo V* 


Consistent 


Coma Ber. 


Bootes I & II 
CVn I & II 
Hercules 
Leo V & Leo T 
Ursa Major I & II 



Table 3: Table of Milky Way ultra-faint dwarfs classified by their distance from our galaxy (columns) and 
whether or not they are consistent with our predictions for the fossils of the first galaxies (row). Note the 
correlation between distance and consistency (*) Pisces II and Leo V are both on the lower end of radii 
expected for fossils, as such they are marked as "inconsistent," but are not as far from predictions as the 
"inconsistent" ultra-faints within 50 kpc. 
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Fi g. 3. — The galactoc entric radial distribution of the fossils excluding the detectable dwarfs as determined 



by IWalsh et al.1 (120091) . Note that the bins have shifted down one order of magnitude in luminosity since 
there are no undetected fossils with Ly > 1O 5 L and we have included the distribution for the lowest 
luminosity fossils with Ly < 1O 3 L . 
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Fig. 4. — The galactocentric radial distribution of the fraction of luminous subhalos which are true fossils 
around MW.2 and MW.3 from Run D. They are divided in the luminosity bins from Figures Q] and |2] for 
which there is a fossil population {Ly < W 6 L@). 
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Fig. 5. — (Left). Cumulative primordial luminosity function of MW.3 from Run D with the observed lu- 
minosity function of Milky Way satellites. We have used a Vfm = 20 km s _1 to determine whether a 
simulated halo is a non-fossil or true fossil. We show the luminosity functions for four distance bins, 
50 kpc< R < 100 kpc (upper left), 100 kpc< R < 200 kpc (upper right), 200 kpc< R < 500 kpc (lower 
left), and 500 kpc< R < 1 Mpc (lower right). In all distance bins the relevant populations are noted as 
follows. The solid black curves show the total cumulative luminosity function from our simulations, the red 
solid lines shows the same for only the star forming halos (non-fossils, including polluted fossils). We show 
the true fossil population with the blue dashed curve. The total observed population is shown as magenta 
triangles with the ult r a-fain t dwarf distribution corrected for sky coverage of the SDSS. The detection limits 
given in lWalsh et al.l (|2009) are shown as vertical black, dashed lines. Long dashed for the luminosity limit 
for the outer radii and shorter dashes for the inner radii in a given bin. (Right). Same a the figure on the left 
but for Vfnt = 30 km s _1 . 
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Fig. 6. — Cumulative primordial luminosity function for MW.3 in Run D and the total observed population 
with only the fossils plotted. The simulated fossils are shown as the blue dashed line and the observed fossils 
as the magenta triangles. The fossil criteria for the observed satellites is the same as in Figure [2] We have 
only shown subhalos around MW.3 which would be detectable by SDSS. 
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Fig. 7. — Cumulative primordial luminosity function f or MW.3 in Run D and observations. Only simulated 
dwarfs which have XV above the lKoposov et all (120081) limit are shown. 
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Fig. 8. — The luminosity (top) and mass (bottom) functions for pre-reionization halos within R v i r of MW.3 
which are not part of a bound subhalo at z = 0. We only include those unbound luminous pre-reionization 
halos between 20 kpc and 50 kpc, where all but two of the known tidal ultra-faints are located. The horizontal 
dotted lines show the approximate number of stripped pre-reionization halos required to reproduce the inner 
ultra-faint population, ~ 30. 
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Fig. 9. — Cumulative luminosity functions of MW.3 from Run D (colored curves) with the total observed 
population (magenta triangles). All the symbols and lines mean the same as in Figure [5] Here we increase 
the stellar mass-to-light ratio by a factor of 10 to 50 M@/Lq. 
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Fig. 10. — Cumulative luminosity functions of MW.3 from Run D (colored curves) with the total observed 
population (magenta triangles). All the symbols and lines mean the same as in Figure [5] Here we increase 
the stellar mass-to-light ratio by a factor of 100 to 500 M©/ 1 L & . 
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Fig. 11. — Cumulative luminosity functions of MW.3 from Run D (colored curves) with the total observed 
population (magenta triangles). All the symbols and lines mean the same as in Figure [5] Here we have 
completely suppressed star formation in any halo outside the highest density regions. This is the most 
extreme case of lower H% formation, and therefore lower star formation in the voids. 
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Fig. 12. — (Left). Cumulative luminosity function of MW.3 from Run D (colored curves) and the total 
observed population (magenta triangles). All symbols and lines mean the same as in Figure [5] We have 
applied an f cr u = 1%. (Right). Same as left panel but for f cr i t = 0.1%. 
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Fig. 13. — (Left) : Radial distribution for our simulated MW.2 and MW.3 with a f cr it = 1% (black solid 
curves) and for the fossil Milky Way satellites (red triangles). Note, that while we are still able to reproduce 
the radial distribution of the lowest luminosity bin, we are no longer able to match the fossil population for 
Ly > 10 5 L . (Right) : Radial distribution for f cr u = 0.1%. Note that for the more extreme suppression, 
we have lost the fossils population with Ly > 10 4 L 
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Fig. 14. — {Left). Histogram of the fraction of luminous true fossils with a given number of luminous pre- 
reionization halos, Ni um . Ni um is a proxy for the number of significant mergers the system has undergone. 
(Right). Histogram of the fraction of non-fossils with a given Ni um . Note, that unlike the Ni um histograms 
for the true fossils and polluted fossils the peak is not at iVj um = 0, but shifted to iVj um ~ 5. Note also, that 
the vertical scale is 0.1 instead of 1.0. 
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Fig. 15. — (Left) '■ Histogram of the fraction of non-fossils with a given V-band surface brightness, £y, 
for the primordial population. We assume that the primordial stars have been puffed up by interactions 
until they fill the full extent of the dark matter halo, giving us £y = Ly x R^ ax , where R m ax is the 
radius at which v(r) = v max . The dashed vertical line is the surface brightness limit of the SDSS from 



Koposov et al.l (|2008T) . Only ~ 1% of the non-fossils are to the right of the dashed line, with expanded 



primordial populations detectable by the SDSS. (Right) : Histogram of the fraction of non-fossils with a 
given V-band surface brightness, £y, for the primordial population. Here we assume that the ghost halo is 
more concentrated, extending only to 0.25 x R max . 
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Fig. 16. — Left : Cumulative luminosity functions of MW.3 from Run D (colored curves) with the total 
observed population (magenta triangles). All the symbols and lines mean the same as in Figure [5] In 
this figure we assume that all the non-fossils are below the detection limits or have lost their entire stellar 
populations due to a combination of heating due to major mergers and tidal interactions. Right : Cumulative 
luminosity functions of MW.3 from Run D (colored curves) with observations (magenta triangles). All the 
symbols and lines mean the same as in Figure [5] Unlike Figure [T6l here we allow the non-fossils to retain 
0.1% of their stellar populations. 



